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We review Floquet formalism of quantum electron pumps. In the Floquet formalism the quantum 
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pump is regarded as a time dependent scattering system, which allows us to go beyond the adiabatic 
limit. It can be shown that the well-known adiabatic formula given by Brouwer can be derived from 
the adiabatic limit of Floquet formalism. We compare various physical properties of the quantum 
O I' pump both in the adiabatic and in the non-adiabatic regime using the Floquet theory. 
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C/3 \ PACS numbers: 73.23.-b, 72.10.Bg, 73.50.Pz 

I. INTRODUCTION 

A quantum pump is a device that generates a dc current at zero bias potential through cyclic change of system 
parameters. The most direct way to create a dc current was originally proposed by Thouless who considered a 
system subjected to a traveling wave potential. This can be realized for example with the help of surface acoustic 
waves 01 ■ Another possibility is to utilize quantum dots. In closed systems operating in the Coulomb blockade 
regime, integer number of electric charge can be transferred by sequential changes of barriers like a turnstile 
whereas in open systems the electron pumping can be driven by adiabatic shape change in the confining potential 
or other parameters which affect the interference pattern of the coherent electrons in the device. In contrast to 
the more familiar dissipative rectification of ac current, the charge transferred in each cycle of adiabatic pumping is 
independent of the period T . After a cycle of the adiabatic shape change we return to the initial configuration, but 
the wavefunction may have its phase changed from the initial wavefunction. This is the so-called geometric or Berry's 
phase 4|. This additional phase is equivalent to some charges that pass through the quantum dot, namely, pumped 
charge |3 . 

i 1 A basic idea of adiabatic quantum pump was theoretically proposed by Spivak et al. 6] in terms of photovoltaic 
effect. It was found that the current linearly proportional to the frequency exists when voltages on the gates U\(t) 
and U2 {t) are arbitrary periodical functions with the same period and different initial phases, and the area enclosed 
by the path in configuration space {U\{i), U2{t)} is not zero. Zhou et al. demonstrated that there exists adiabatic 
quantum pumping different from the usual photovoltaic effect. Both the amplitude and the sign of the charge 
transferred through a sample per period are random sample spcecific (i.e. not quantized) determined by the quantum 
, interference. Switkes et al. experimentally have realized such an adiabatic quantum pumping through an open 
• quantum dot under the shape deformation controlled by two ac gate voltages At the same time Brouwer [fj 
proposed simple and intuitive approach to deal with the quantum pump using adiabatic scattering matrix, which can 
explain many aspects of the experimental results perfomed by Switkes et al. . Since then, adiabatic quantum pumps 
-t— > ' have attracted considerable interest from many researchers [l0| . 

In other sense, the quantum pump is a time dependent system driven by (at least) two different time periodic 
perturbations with the same angular frequency and initial phase difference 4>. One can deal with this problem using 
not only adiabatic approximation exploited by Zhou et al. and by Brouwer but also the so-called Floquet approach 
[Tl|. Floquet method allows us to solve the time dependent (mostly periodic) problem precisely. In physical ponit of 
view, an oscillating potential can transfer an incoming electron of energy E to Floquet side bands at E± nhio, where 
n is an integer and uj is the angular frequency of the oscillation. A scattering matrix for a time dependent system can 
be constructed from the interplay of these sidebands 0, ^| . The Floquet approach can bring us beyond adiabatic 
regime of quantum pumps. In this review we will present a breif review on the floquet formalism of quantum pump 
developed recently. 
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II. FLOQUET FORMALISM FOR TRANSPORT PROBLEM 

Moskalets and Biittiker [l4|, and Kim have developed the Floquet scattering theory of quantum pump inde- 
pendently First, we breifly introduce the derivation of pumped current under the Floquet theory. 
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A. Directed charge currents 



Consider the one-dimensional time-dependent Schrodinger equation ih(d/dt)tp = H{t)tp for a non-interacting elec- 
tron with mass fi and H (t) — — h 2 W 2 /2fi + U(x, t), where U(x, t + T) = U(x, t) and U(x, t) = at x — > ±00. Due to 
the periodicity in time, a solution can be written as 



* t {x,t)=e- ut ' h 



(1) 



where e is the Floquet energy which takes a continuous value in the interval [0, hw). Since the potential is zero at 
x —>■ ±00, Xn(x) is given by the following form 



Xn{x) 



A n e ik " x + B n e- ik " x , x -> -00 
C e ik n x + n e -ik n x x ^ +QOi 



(2) 



where k n = y/2fi(e + nhuj)/H. By wave matching, the incoming and the outgoing waves can be connected by matrix 
M 



= M 



(3) 



If we keep only the propagating modes (k n is real), we can obtain the uni tary Floquet scattering matrix [see the 
argument below Eq. (|A8() ]. which can be expressed in the following form [lallallq ] 
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where r a p and t a p are the reflection and the transmission amplitudes respectively, for modes incident from the left 
with an Floquet energy e which take continuous values in the interval [0, hui); r' a p and t' a p are similar quantities for 
modes incident from the right. 

To get a current let us consider the scattering states. For an energy E = H 2 k 2 /2/j, (k > 0) of the incoming particle 
the scattering states as a solution of the Schrodinger equation can be defined as 



X E (x,t) = 



gikx — iEt/h 

Efi n >0 ^E n E 
p —ikx — iEt/h 
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+00, 

+00, 
—00, 



(5) 
(6) 



where E n = E + nHui, k n — ^/2fiE n /h, and the normalization is ignored. Here we have reflection and transmission 
coefficients r E E and t E E , which can be obtained from the unitary Floquet scattering matrices S obtained above. 
\Hi El ■ The above transmission and reflection coefficients are related to the matrix elements of S in terms of 

+ 



t 



E„,E 



y/kp/k a t a p> etc. with E — e + fihw and a = n + (3. 



We derive the current u sing these scattering states. The time dependent electron field operator can be obtained in 
the following form [13, 13 Ei 



O-trE 



(7) 



where a a E and v(E) is an annihilation operator for electrons in the scattering states Xe( x ^) ano - the velocity, 
respectively. Even though the scattering states do not form othogonal bases we need only the completeness to be sure 
that the expansion of Eq. (JJJ is valid. Using this field operator the current operator is also expressed as 



J(x,t) = (ie/2m)* + (x,t)V*(x,t) + H.c. 



/ dEdE' Im( X °E,*Vx E ) a ° Eaa ' E . 
J V E hy/v(E)v{E') 



(8) 
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The quantum mechanical (or thermal) average of the current operator becomes 

cr (7(7 i? re >0,n7^0 V V / V 

We evaluate Eq. (jSJ) taking x — ► oo and averaging over space and time. One can then obtain the pumped current as 
following 

7 = \ E / - T E nE fR(E)] , (10) 

E n >0 J 

where we exploit X)_E„(^«/^)l r B _el 2 = 1 — X^^n/^OI^fil 2 an( l ( a aE a <?E) = fa{E). The second term of the 
righthand side in Eq. jjjjj vanishes due to the unitarity of the scattering matrix j^. Here, T E E denotes (k n /k)\t E E \ 2 . 
The hnally obtained current expression simply means that the pumped current corresponds to the difference 
beween the currents going from the right to the left and from the left to the right. 



B. Adiabatic limit 

Moskalets and Biittiker 0] have shown that the adiabatic limit of the Floquet formalism is equivalent to the 
adiabatic scattering matrix approach given by Brouwer. Here we summarize it. 

The adiabatic condition in the quantum pump implies that any time scale of the problem considered must be much 
smaller than the period of the oscillation of an external pumping 9] . We can then define the instantaneous scattering 
matrix with time dependent parameters, namely X n (t), 

S ad (E,t) = S ad {E,{X n {t)}) = ( Y ) (11) 

\ l ad r ad J 

Due to the time periodicity of X„'s, using a Fourier transform one can obtain the amplitudes of side bands for particles 
traversing the adiabatically oscillating scatterer with incident energy E as following 

S ad (E, {X n (t)}) = S ad , n {E)e- iw \ (12) 

n 

where 

1 f T 

S ad , n (E) = - dte m ^S ad (E,{X n (t)}). (13) 
1 Jo 

Thus we can construct the adiabatic Floquet scattering matrix as follows 

S{E n , E) » S{E, £_„) = S ad ,n(E), (14) 

where E n = E + nfkv. 

Now we go back to the current expression Eq. IjlUI) in the Floquet formalism. Assuming Jl(E) — fn(E) = f(E), 
Eq. (fTUf> can be written as 

1 = I E f dE [\ s kE\ 2 m - \s EB _ff(E- n )] , (is) 

where we make the shift E — > E — nhui for the second term in the bracket. Here one can see T E , E = \S^, E \ 2 since 
tnm(t'nm) = ^e'e(^E'e) with E' = e + Tifiuj and E = e + mfajJ. Using the adibatic approximation Eq. (|14Jl and 
lu << 1, we obtain 

™ f^df(E) 



Iad = Y dE dE E n\S ad AE)\ 2 - (16) 

77 E„>0 



Using Eq. (|12|) . Eq. (|16|l can be rewritten as follows: 



T 



'4tt 2 ,/ n J n ' dE " \ at ad 

This is Brouwer's formula. 
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FIG. 1: Even though the system has time reversal symmetry (S = S T ), the transmission probabilities for the opposite directions 
can differ from each other. Note that in Floquet formalism the transmission probability is obtained from taking summation of 
one column. 



C. Quantum pumping and broken symmetry 

In the adiabatic limit the pumped current occurs when the area enclosed by the path in the parameter space 
{Vi(i), V2(f)} is not equal to zero. It implies that the system breaks the time reversal symmetry. The time reversability 
of the time dependent system is subtle. We define that the system is time reversible if the Hamiltonian satisfies 
H(t) = H(—t) after appropriate time translation H(t + a) with a certain a. For simplicity let us assume the potential 
V(r, t) is given by 

V(r,t,4>) = Vi(r)cos(wf-0/2) + V 2 (r) cos(ojt + 0/2), (18) 

where V\ and V2 represent the spatial dependence of two time dependent perturbations. <p is t ne initial phase difference 
of two time-dependent perturbations. In case (f> = one finds that 

V(r, t, 4>) = [Vi (r) + V 2 (r)] cos tot, (19) 

which means effectively there is only one time dependent perturbation, so that no pumped current exists. Even 
though the cosine function is replaced by sine, by adding 7r/2 to the time t, one can make the system retain time 
reversal symmetry. It is emphasized that in the adiabatic limit the time reversal symmetry should be broken in oder 
to obtain the pumped current. 

In the non-adiabatic case, however, the condition that the pumped current appears is quite different from that in 
the adiabatic limit. From Eq. (|10|) one can see the pumped current exists when T E E ^ T E E is fulfilled. We assume 

that the system has the time reversal symmetry. The time reversed process of the transition Tg E is given by T E E 
since not only one should change the directions of the momenta but also replace the emission by the absorption, and 
vice versa. It means that for the system with the time reversal symmetry T E E = T E E holds. This is nothing but 

S T — S, where T represents the transpose of a matrix. However, such an equality does not guarantee Tg E = T E ^ E . 
In other words, the pumped current in the non-adiabatic case can still exist even when the time reversal symmetry 
holds. In the adiabatic limit (namely, ui — ► 0), T E E approaches T EE ■ Then the time reversal symmetry, i.e. 
Te e„ = T E E , leads us to T E E w T E E , i.e. no pumped current. This is the reason why the pumped current 
vanishes in the adiabatic limit when the time reversal symmetry is preserved. 

To get more physical intuition for the relation between the pumped current and the broken symmetry let us 
reconsider the current expression given in Eq. l|lUfl • The pumped current / is given by the difference of two currents 
having the opposite directions. The transmission probability to the right for an electron with the incoming energy E 
is ~^2 E > qT e e (= T E ). Comparing it with Eq. Q one can find that T E correponds to the summation of the mth 
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column of a part of the Floquet scattering matrix when E = e + mfku (m is an integer). More precisely ^2 E >Q Tg E — 
^^l t nm . It is worth mentioning that in usual multi-channel scattering problems with time-independent potentials 
the transmission probability is given by the summation of a block of a scattering matrix, e.g. ^2 mn t nm , where m 
and n represents the number of modes or scattering channels of the incoming and the outgoing lead, respectively. 
Even though the system considered has time reversal symmetry, i.e. S = S T or equivalently t nm — t' mn , it is not 
guaranteed that J2n *nm = Y, n ( see Fi § wnile & is a l wa ys true that J2 nm ^nm = J2, lm t'nm once t mn = t' mn 
is fulfilled, which is related to the reciprocity. The physical meaning of the summation of one column in the case with 
the oscillating potential to obtain T E is the following. The incoming electron has a definite incident energy E since 
the energy is well defined far from the scatterer, whereas the outgoing electron can have various energies to be taken 
into account since an oscillating potential can transfer the electron with energy E to side bands at E ± hui. 



D. Wigner delay time 

The adiabatic approximation implies that any time scale of the problem considered, especially the electron dwell 
time in a quantum pump must be much smaller than the period of the oscillation of a external pumping T p [jj. 
Using the Floquet formalism we can calculate the Wigner delay time tw , which is the interaction time of the incident 
electron with the scattering potential 0, |2^, . In this sense Tw corresponds to the electron dwell time in the 
quantum pump. To obtain the Wigner delay time we use the eigenvalues of the Floquet scattering matrix 5. Due to 
the unitarity of S all the eigenvalues lie on the unit circle and can be written in the form exp(i8 a ). The Wigner delay 
time is defined by 

T w = hJ2^\(k n \e a )\ 2 , (20) 

a 

where the eigenstate corresponding to the eigenvalue 9 a and an input propagating state (or channel) with momentum 
k n are denoted by \9 a ) and \k n ) respectively 24]. It is worth noting that the Wigner delay time is a function of the 
energy of the incident particle E (= Epi + nhu), and \k n ) and \9 a ) are determined by n and Epi respectively. If 
(k 7l \9 a ) is ignored in Eq. 1)20(1 the Wigner delay time tw becomes trivial, i.e. tw(E + nftuS) = Tw{E). 



III. PAULI BLOCKING FACTORS 



In mesoscopic systems it is well known that the current through the scatterer can be obtained from 

1=1 [ dEdE' [T+(E',E)f L (E)-T-(E,E')f R (E / )] , (21) 

where T + (E',E) represents the transmission probability for scattering states incident from the left at energy E 
and emerging to the right at E' , and T^{E 1 E') is defined in a similar manner for the reverse direction, {fii) 
is the Fermi-Dirac distribution in the left (right) reservoir. There has been some debate on using this formula 
mmmmiMIMIIllii since E< 1- ® dose not contain the so-called Pauli blocking factors. The fermionic 
nature of the electrons is taken care of in an ad hoc way by factors 1 — / to suppress scattering into occupied states, 
so that the current is given by 

7= e - J dEdE' {T+(E',E)f L (E) [1 - f R {E')\ - T~(E, E')f R (E r ) [1 - f L (E)]} . (22) 

Usually these two expressions Eqs. (|21|) and 12 I'D give the same results since the difference between them, [T + (E' , E) — 
T-(E,E')]f L (E)f R (E'), vanishes when T+(E\E) = T~(E,E'), i.e., the micro-reversibility holds. The question can 
arise, however, if the system lack of this micro-reversal symmetry is considered. One of the relevant example is a 
quantum pump. 

More than two periodically oscillating perturbations with a phase difference (7^ nn, n is an integer) break the time 
reversal symmetry 32] as discussed in Sec. II. C, and consequently T + (E',E) ^ T~( E, E'). Therefore, the currents 
obtained from Eq. (|21|l and 11221) are different from each other in quantum pumps [29t |32( . The question immediately 
arises: which one is correct in quantum pumps? It is noted that this problem still exists even in the adiabatic limit. 
We would like to make a conclusion first. Even though the time dependent scatterer like the case of quantum pumping 
not only cause inelastic scatterings but also break the mico-reversibility, the Pauli blocking factor is unnecessary when 
the scattering process through the scatterer is coherent. 
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The existence of the Pauli blocking factors is intimately related to the "scattering states" [25|. If we fill up the 
energy eigenstates with the electrons in both reservoirs independently and then transfer the electrons from one to 
the other reservoir, the Pauli blocking factors cannot be avoided. If the transport is coherent across the scatterer, 
however, one can define a single wavefunction extending from one reservoir to the other (more precisely reflected 
and transmitted waves in every connected reservoirs) and then fill up these scattering states. In this consideration 
the concept of transferring the electron from one to the other reservoir is automatically eliminated, so is the Pauli 
blocking. In Sec. II. A using the scattering states given in Eq. JSJ and © we have already shown that the current 
experssion ()10|l does not contain Pauli blocking factor. 

Recently, Wagner has raised the question on the Pauli blocking factor in quantum pumps, and concluded that a 
marked difference is predicted for the temperature dependence of the pumped current for the case with and without 
the Pauli blocking factor |32j . The pumped current can be generated from two contributions: classical and quantum 
origin. Most interest has been focused on the quantum mechanical charge pumping in chaotic quantum dots, where 
the pumped current has zero average over dot configurations and shows the mesoscopic fluctuation [tj. The 

1/2 

important measure quantifying the currents pumped quantum mechanically is a mean square average (/ 2 ) . Since 
the mesoscopic fluctuation of the pumped current results from the electron interfence and is generically quantum 
nature, the scattering process in physically meaningful quantum pumps should be coherent. It was observed in 
experiment || that as the temperature increases the mean square average decreases due to the temperature-dependent 
decoherence. The Pauli blocking factor is then unnecessary from the beg inning, which is also consistent with the results 
obtained from non-equilibrium Green's function formalism in Ref. [29j. 



IV. MAGNETIC FIELD INVERSION SYMMETRY 

There has been the puzzle unsolved in the experiment performed by Switkes et al. @, which is the magnetic 
field inversion symmetry (MIS) in pumped currents. It was suggested theoretically Q that the pumped current / is 
invariant upon magnetic field reversal 

1(B) = I(-B). (23) 

The subsequent experiment on open quantum dots appears to be in a good agreement with Eq. (|23|) 8]. It became 
immediately clear, however, that such symmetry is not valid in theory |33|. In addition, the measured dc voltage 
is 20 times larger than the theoretical estimation. Furthcrmocre, the voltage fluctuation decreases upon applying 
the static magnetic field breaking time-reversal symmetry, which seems to be the weak localization-like phenomenon 
of disordered systems, while such an effect cannot be predicted in the theory. Brouwer proposed that the rectified 
displacement current could account for all these discrepancies |34j . which implies the dc current observed in the 
experiment is not attributed to the adiabatic quantum pumping. Recently Marcus grou p has reported that the origin 
of the dc current is mesoscopic rectification as well as adiabtic quantum pumping [35j. Strictly speaking, a reliable 
adiabatic quantum pumping has not been realized yet. 

Even though in general quantum pumps do not fulfill the MIS, additional discrete symmetries of quantum dots 
can lead to such MIS. The effect of discrete symmetries on the MIS has been studied by Aleiner, Altshuler and 
Kamenev |36| in the adiabatic limit. They found the reflection symmetries give rise to relations 1(B) = I(—B) or 
1(B) = —I(—B) depending on the orientation of the reflection axis. In the presense of inversion 1(B) = 0. Note 
that the symmetry considered in Ref. [3^ | should be kept intact in the pumping cycle, and the results obtained are 
available only in the adiabatic regime. Using the Floquet approach we can also investigate the MIS of the pumped 
current for three discrete symmetries, namely LR (left-right), UD (up-down), and IV (inversion), as shown in Fig. 
both in the adiabatic and in the non- adiabatic cases. These symmetries are only applied to the time independent part 
of the scattering potential, which is V n (r) in Eq. 1251) , with an arbitrary cf>, so that they are not necessarily kept intact 
during the pumping cycle |37j . 



A. Floquet approach 

Let us consider the two-dimensional quantum dot driven by time-periodic perturbations under a static magnetic 
field, which is attached to two leads (see Fig. whose Hamiltonian is given by 



..dip 



(ihV + eAf 
2^ 



+ V(r,t) 



1>, 



(24) 




FIG. 2: The schematic diagram for three discrete symmetries of quantum pumps, i.e. (a) LR, (b) UD and (c) IV symmetry. 



where A is the vector potential. The potential V(r, t) is given by 

V(r,t, 0) = Vb(r) + Vi(r) cos(^ - 0/2) + V 2 (r) cos(u;f + 0/2), (25) 

where Vo represents the confined potential of a quantum dot, and V\ and Vi are the spatial dependence of two time 
dependent perturbation. is the initial phase difference of two time-dependent perturbations, and r denotes (x, y). 
Note that V(r, — t, 0) = V(r 7 t, —(/)). When — nir without magnetic field, the quantum pump is microscopically 
reversible (or time reversible), i.e. = where represents the transmission amplitude from the Floquet side 
band m of the channel a in the left lead to the side band n of the channel (3 in the right lead, and t'£& is a similar 
quantity for the opposite direction. In a quantum pump the two time dependent perturbations with a finite nit) 
as well as the magnetic field break the time reversal symmetry as already discussed in Sec. II. C 

If we take the complex conjugate of Eq. I|24l) and at the same time, reverse the vector potential (A — > —A) and 
time (t —>■ —t), we obtain 



6V* 
~dt 



^ + eA) \v(r,-t) 

2[A 



ip*. (26) 



Compared with Eq. (|24|) one can see that i/j-b,$(y, t) — tpg _^( r , — t). The complex conjugation combined with the 
time inversion of the solution in Eqs. Q and © simply corresponds to the reversal of the momentum direction of the 
incoming or the outgoing plane waves and simultaneously the replacement of absorption by emission, and vice versa. 
This simply implies that S — > S T (S T is the transpose of S), i.e. — > S^ n . For completion of the magnetic field 
inversion operation, 4> — > —<fi also has to be considered. 

The total transmission coefficient to the right reservoir after the magnetic field inversion is given by 

oo 

T -bA E ) = EEI*nmfe-S,0)| 2 (27) 

otB n=0 
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FIG. 3: Schematic diagrams for describing the configuration of side bands with given channels a and (5 to calculate the total 
transmission coefficients through the time-periodic perturbation V(t) for (a) T^AE) and (b) the case given in Eq. 1281 , (c) 
The setup (a) (the solid arrows) superimposed by the setup (b) reversed by using IV symmetry given in Eq. do 1 1 (the dashed 
arrows) . 



= EEIC^,5,-0)| 2 , ( 28 ) 

a/3 n=0 

where E = e + mhw. It must be noted that Eq. I|28|l is not equivalent to Tg _,. Figures Eta) and (b) represent the 
usual setup of side bands for calculating the total transmission coefficients, e.g. AE), and the setup for calculating 
Eq. 128(1. respectively. In general, they are different from each other, so are 1(B) and I(—B). The MIS is not valid in 
Floquet formalism, as well. 

LR symmetry. If we assume LR symmetry [Vb(x, y) — Vo(—x, y) and Vi(x, y) — V2(— x, y)\, the following relations 
are obtained [see Fig. Efa)] 

t^(e,B,<f>) = t'^(e,-B, -</>), 

r^{e,B,4>) = r'^(e-B,-<f>). (29) 
These are expected from the symmetries of classical trajectories (3?|. Eq. (|27f) can then be rewritten as 

oo 
a/3 n=0 

= T B ^(E). (30) 

From Eq. QlOfl one reaches I(B, <f>) = —I(—B, —<f>). Since such relation results from the Floquet approach, it is available 
for both the adiabatic and the non-adiabatic cases. At zero magnetic field one finds I((j>) — —/(—</>), consequently 
I(4> — 0) = even for the non-adiabatic case. 

UD symmetry. For UD symmetry [Vo(x,y) = Vo(x,—y) and V\[x,y) — V^ix,— y)], from the similar procedure 
used in Eq. (|29|l and (|3L)[I one finds I(B,cf>) = I(— B,— 4>). Like the LR symmetry, such symmetry of the quantum 
pump is correct for both the adiabatic and the non-adiabatic cases. At zero magnetic field one finds 1(4') — I{~ 
which can be understood from the fact that the exchange of V\ and V%, i.e. <f> ~ * ~4>i does not make any difference in 
Fig. lib). 

IV symmetry. If one assume IV symmetry [Vo(x,y) = Vo(—x,—y) and V\(x,y) — Va(— x, — y)], one obtains [see 
Fig. He)] 

tfa l (e,B,<f>)=ltf£{e,B, -<j>), 

r^(e,B,<fi)=r'P°(e,B,-<j>). (31) 

Eq. is reexpressed as 

oo 

T -^( £ )=EEl i lfcB»| 2 , (32) 

a/3 n=0 
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which is different from Tg ^(E) since the summation is taken for n as shown in Fig. OJc). In general, there is no 
relevant MIS for dots with IV symmetry. When tj> = Q, however, Eq. (|31fl leads to Tg(E) = Tg(E), consequently 
/ 0. 

B. Adiabatic limit 

We have already investigated the relation between the Floquet and the adiabatic approach in Sec. II. B. Once the 
adiabatic scattering matrix is found, by using Brouwer's formula the pumped charge Q per one cycle can be 
expressed as 

Q{B) = g[i ad (B),r' ad (B)} = f dt{f[t ad (B)} + f[r' ad (B)]h (33) 

Jo 

where 

k=l a/3 

Then, the pumped current I is Q/T. Below, we will use Eq. (|33[) for the current. We introduce I r and which 
are the pumped current to the right and the left, respectively, so that I r = — Ij. I r is nothing but / in the previous 
notation. 

LR symmetry. The relation Q29J1 implies i(r)^ d (E;B,<f>) — i'(f')^ dn (E;—B,—(f)) in the adiabatic limit, which, 
by using Eq. l(T2|) . immediately leads to i(f) ad (E, t; B, (f>) — t'(f') a d(E, t; —B, —(f)). From Eq. one finds 

I r (-B,<p) = g[i{-B,4,),f'{-B,4>)] 

= g[i'(B,-ct>),r(B,-ct>)} (35) 
= Ii(B,-<f>) = -I r (B,-cl>). 

Note that if the relation I(<fi) — — 1(— (j>), which is available only in the adiabatic limit, is considered, one reaches 
simpler form 1(B) = I(—B). 

UD symmetry. Like the previous subsection one can obtain t(r)^ d n (E; B , <p) — i(r)^°l n (E;—B,—<j)) in the 
adiabatic limit, which, by using Eq. (H21, leads to i[f) ad {E, t; B, 4>) = i(f) ad (E, t; -B, -<j>). From Eq. O one obtains 

I r {-B,cj>) = g[t{-B,(l>),r'{-B,4>)} 

= g[t(B, -<(>), f'(B,-cf>)] (36) 
= Ir(B,-<P). 

Note also that if the relation 7(0) = — /(— <j>) is considered, one reaches simpler form 1(B) — —I(—B). 

IV symmetry. The relation (|3*T1) implies i(r)^ n (E: B,<j)) = i' (f')^ d n (E; B, —<f>) in the adiabatic limit, which, by 

using Eq. (fT2|> . leads to i(f) ad (E 7 1; B, cj)) = i'(f') ad (E, t; B, —0). Considering the magnetic field inversion operation, 
one can obtain 

t ad (E,t;B,4>) = F ad (E,f t -B,<j>), 

r ad (E,t;B,<f>) = fZ(B,t;-B,cj>). (37) 

From Eq. I|33|l one finds 

I r (-B,4>) = g[t{-B,cl>),f'(-B,4,)} 

= g[F(B,4>) 1 r T (B,4>)]- (38) 

There is no relevant symmetry in the adiabatic limit, again. 

Eq. (|14fl implies that t m +n,m ~ tmm-n m the adiabatic limit. After summing up all the side band contributions, 
the two configurations shown in Fig. [21(c) (the solid and the dashed arrows) give the approximately equivalent total 
transmission coefficients, i.e. Tg(E) w T^ B (E). By applying similar procedure to the reflection coefficient one can 

also find R^(E) ps RT b (E) (make sure that the arrow is now reversed). It is also satisfied that Tg(E) sa TZ B (E) 
since the unitarity of Floquet scattering matrix requires T + = 1 . From these two relations the pumped current 
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TABLE I: The symmetries of the pumped currents for the magnetic field inversion. TRS denotes time reversal symmetry by 
two time dependent perturbations, i.e. TRS for (j> = nir and no TRS for <f> 7^ nir. x represents that there is no relevant 
symmetry. 



adiabatic 






non-adiabatic 




TRS no TRS 




TRS 


no 


TRS 


LR I = 1(B) = I(- 


B) 


1(B) = -I( 


-B) I(B,<f>) = 


-I(-B,-<t>) 


UD / = 1(B) = -I(- 


-B) 


1(B) = /(- 


-B) I(B,f) = 


■ I(-B,-4>) 


IV 1 = 7^0 




1 = 




X 



is shown to be vanishingly small, 1(B) = J dE[T+(E) - Tg(E)] » since T+(E) » T+ B (E) » Tg(E). In Table I, 
we summarize the MIS of the pumped currents with various situations for the three symmetries. 

One possibility to realize the experimentally observed MIS, i.e. 1(B) = I(—B), is that the quantum dot would obey 
LR symmetry. Then, the MIS is recovered only in the adiabatic limit, while it will be broken in non-adiabatic case. 
In general, however, it is difficult that a chaotic quantum dot possesses any discrete symmetry like the LR symmetry. 



V. SIMPLE EXAMPLE: TWO OSCILLATING DELTA-FUNCTION BARRIERS 

The minimal model that shows quantum pumping effect is one-dimensional (ID) two harmonically oscillating delta- 
functions barriers. Wei et al. have applied Brouwer's adiabatic approach to this model to calculate the pumped current 
|39| . Using this simple model one can investigate many important properties of quantum pumps. 



A. Brouwer's adiabatic approach 



We consider ID two harmonically oscillating (5-function barriers with the strengthcs X\ = Y\ + Aicoswr and 
X 2 = V 2 + A2 cos(ojt + (f>) respectively, separated by a distance d. This is a simplified model of the experiment by 
Switkcs et al., but possesses many important characteristics and can be easily handled. Wei et al. studied parametric 
charge pumping aided by quantum resonance using this model and found that the pumped current has large values 
near a resonance level |39j. Due to the double barrier geometry, resonant tunneling also plays an important role in 
charge pumping. This is not so surprising because the adiabatically pumped charge is proportional to the variation 
of the scattering matrix to the external parameters. Usually the elements of a scattering matrix are very senstive to 
the variation of the external parameters near the resonance. 

The 2x2 scattering matrix of the double barriers with the strengthes X\ and X 2 is given by 

*- (:")■ <») 

where r and t are the reflection and the transmission amplitudes respectively, for modes incident from the left; r' and 
t' are similar quantities for modes incident from the right. The charge emitted per cycle to the right is obtained from 

f Tp , ( dm dX 1 dm dX 2 \ 







dXi dr dXi dr 



where 



dn% T { dt iSf dr' lJf 



dx m - lia {dxZ r + dxZ r J (41) 

(m = 1,2), and T p (— 2tt/uj) is the period of the pumping. One can show that the emitted charge to the left Q2 is 
equal to —Q\. Equation 140(1 can be rewritten in the following form by using Green's theorem 

Qi = e / dX 1 dX 2 U(X 1 ,X 2 ), (42) 

J A 

where H(Xi,X 2 ) = d(dn\/ dX 2 ) / dXi — d(drii/dX 1 )/dX 2 - The pumped current is easily obtained from Zj = Qi/T p . 
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FIG. 4: The Wigner delay times (a) with the same condition as (c), and (b) as (d). The pumped current 7i calculated by using 
Brouwer's approach (solid curves) and the Floquet approach (dotted curves) with <j> = n/2 for (c) A = 22.5 meV-nm and (e) 
A = 225 meV-nm with V = 0, and (d) A = 22.5 meV-nm and (f) A = 225 meV-nm with V = 225 meV-nm. In (d) and (f) the 
transmission resonances are denoted by the filled circles (•) obtained from considering static double barriers, whose y values 
are chosen arbitrarily. The attached error bars represent the sizes of the imaginary energy of each resonance. 



We use the parameters d — 50 nm for the distance between the two barriers, the effective mass [i = 0.067m e of an 
electron in GaAs, and T p = 9.09 ps for the period of pumping, which corresponds to fkv = 0.45 meV. Figure 01 shows 
the pumped current as a function of the energy of an incident electron with V\ — Vi — V, Ai = A2 = A, and <j) = tt/2. 
We present two examples; the first is a nearly open case (V=0) in Figs.^Ji and0^, and the second is a closed case 
(V=22h meV-nm) in Figs. and^f. All of them show interesting resonance-like structures (See Ref. 0] for the 
details of the resonance-like behavior). 



B. Floquet approach 

Now we study the Floquet approach of the problem investigated in Sec. II. Using the scattering matrix of a single 
5- function with sinusoidal time dependence (see Appendix for details), we can obtain the total scattering matrix of 
the oscillating double (5-functions in the following form 

-(:::)• 

where r = tl + t/,(I — Qr i? Qr £ ) _1 Qr i? Qt i , and t = t#(I — Qr L Qr^) _1 Qt £ ; r' and t' can also be obtained by 
replacing L by R in r and t respectively. Here Yr(l) an d ^r(l) are the reflection and the transmission matrices 
respectively, for the right (left) delta function with time dependence for modes from the left, and I is an identity 
matrix. Due to the reflection symmetry of each delta function y' r , l s = i\r(z,), and = ^r(l)- During each one-way 
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trip an electron at energy E n — E + nhu> picks up a phase factor exp(ik n d), which is represented by the diagonal 
matrix Q m „ = exp(ik m d)6 mn . Once we have Floquet scattering matrix S, the pumped current to the right is given 
by Eq. jfOJ), i.e. 

I = | J dEdE'[t(E',E)f L (E) - t'{E,E')f R (E')]. (44) 

Without external bias (fi = fn = f) Eq. (|44(l can be rewritten as 

2e f 00 

/=-y o dEf(E)[T + (E)-T-(E)}, (45) 
where T + (E) = >Q T En E . At zero temperature it becomes 

I=\ ' dE[T+(E)-T-(E)}, (46) 

where Ef is the Fermi energy. 

Figure 0] shows that in the open case (V = 0) for small A (Fig.^Ji) the pumped current obtained from the Floquet 
approach is equivalent to that of Brouwer's while for rather larger A (Fig. 0^) they deviate from each other, which is 
not fully understood yet. When V ^ (the closed case), even for small A (Fig.^Jl) they are quantitatively different 
near the resonances. It is answered in the following subsection. 



C. Wigner delay time 



Figures^ arid b show the Wigner delay times by using the same parameters exploited in Fig.^J: and d, respectively. 
In Fig. 2Ji the Wigner delay time decreases as E increases, which can be understood when one takes it into account 
that usually the electron dwell time is short if the energy of an incident electron is large. Note that the sharp spikes for 
low energies are just artifacts coming from the fact that the Floquet scattering matrix numerically obtained diverges 
at E — nfojj. Since the Wigner delay time is small enough, the adiabatic condition is still fulfiled to show very good 
agreement between the two currents in Fig.QJ;. In Fig. 2Ja several peaks of the Wigner delay time is clearly shown. It 
can be explained by thinking that now two static barriers are plugged into the problem, which generates the so-called 
double barrier resonances. The Wigner delay times have large values near the resonances since at the resonances an 
electron can stay in the quantum dot for a long time, which immedeatly implies that the adiabatic condition can be 
broken down. It explains the deviation between the pumped currents near the resonances observed in Fig. 0Ji. 



D. Broken symmetry 

One of the interesting consequences from Eq. 1|46|) is that the pumped current still exists even in the cases <f> = or 
7r when Ai ^ A2 ^3- Since the integration area in parameter space is zero when <fi = or tt, in Brouwer's approach 
the pumped current definitely vanishes. In contrast, even when <f> — or 7r, an asymmetry of the potential can lead to 
the asymmetry of the currents which is nothing but the pumped current in Eq. (|46|l . Figure shows the pumed 
current as a function of ratio of the strength of two barriers A2/A1 with <f> = or tt, and E = 6.005fiu;. Note that the 
pumped current is zero when Ai = A2. The oscillatory behavior is also related to the double barrier resonances. 



E. Pauli blocking factor 

Figure shows the pumped currents obtained from Floquet approach with and without Pauli blocking factor, and 
Brouwer's formula in the same model used above under the adiabatic regime as shown in Fig. ^1. It is clearly seen 
that the current with Pauli blocking factor deviates from that of Brouwer's approach which nearly coincide with the 
current without Pauli blocking factor. It is worth nothing that qualitative behavior of the pumped currents with and 
without Pauli blocking looks quite similar: the pumped current / oc A 2 sin(/) as shown in the insets of Fig. 
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FIG. 5: Pumped current I\ as a function of potential asymmetry A2/A1 with Ai = 22.5 meV-nm and E = 6.005/k<j at <f> — 
(solid curve) and <j> = n (dashed curve), where the currents are normalized to their values at at = n/2. Inset shows the 
magnification of a part of the plot. 
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FIG. 6: The pumped currents obtained from the Floquet approach with (the thick solid curve) and without Pauli blocking (the 
thin solid curve), and Brouwer's formula (the dashed curve) with the same parameters used in Fig.0 Note that Io — e/T = 17.6 
nA. The upper inset: the pumped currents as a function of the strength of the oscillating potential A from the Floquet approach 
with (the filled circles) and without Pauli blocking (the open circles), and Brouwer's formula (the solid curve). The lower inset: 
the pumped currents as a function of a phase difference <j> from the Floquet approach with (the solid curve) and without Pauli 
blocking (the dashed curve). In both insets we take E = 6.005/kj 
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VI. SUMMARY AND DISCUSSION 



In summary, we have reviewed Floquet formalism of quantum electron pumps and shown its various application. 
In the Floquet formalism the quantum pump is regarded as a time dependent scattering system, which allows us 
to go beyond the adiabatic limit. We have derived the expression of the pumped current using scattering states, 
and shown that the adiabatic formula given by Brouwer is successfully recovered in the limit u — > 0. The condition 
for obtaining non-zero pumped current is rather subtle when the non-adiabatic case is considered. In the adiabatic 
case the time reversal symmetry should be broken to have the pumped current, while it is sufficient to generate the 
pumped current to break either the spatial reflection symmetry or the time reversal symmetry in the non-adiabatic 
case. Since usually the quantum pump can break the time reversal symmetry with time dependent perturbations, 
the issue of Pauli blocking originating from the Fermionic nature of electrons can occur. It can be shown, however, 
that the Pauli blocking factor is not neccessary at all. Using the Floquet formalism one can discuss the symmetry 
properties of the pumped current over the inversion of the static magnetic field for various spatial symmetries of 
quantum dots. Application of the Floquet formalism dierctly to realistic problems such as 2D quantum dots is rather 
difficult because one should consider many sidebands for each scattering channel. However, we can get many basic 
important physics of quantum pumps by studying a rather simple model system: a ID oscillating double (S-function 
barriers. 

The adiabatic charge pumping does not have any corresponding classical analog since it has something to do 
with Berry phase, which means this is completely quantum mechanical phenomena. The quantum pump in the 
non-adiabatic regime, however, is somewhat subtle because Berry phase concept dose not work any longer. It has 
been shown that there are several other mechanisms to generate dc current through the non-adiabatic time periodic 
perturbations with zero average, e.g. quantum ratchets [44L l45l |46| . quantum ballistic rectifiers [471 l48j . Hamiltonain 
ratchet , and quantum shuttles [50, • We believe that it must be important to understand the relations and 
the differences among these devices. The Flouet formalism might play an important role to understand them because 
it allows us to go beyond the adiabatic limit of quantum pumps. Finally, we would like to mention that it is also 
important to know what role the decoherence plays in the system with time periodic perturbation. It helps us 
understand the classical-quantum correspondence of dc current generation of quantum pumps, ratchets, rectifiers, 
and shuttles. 

The Floquet formalism has been successively used to describe a lot of physical properties of the quantum pump. 
It is a powerful tool not only for understanding quantum pumps, which is the main topic of this review, but also for 
many kinds of coherent transport problems with time-periodically varying scatterers. 



This work was supported by Research Institute for Basic Sciences, Pusan National University (2004), and Pusan 
National University Research Grant, 2004. 



APPENDIX A: FLOQUET SCATTERING MATRIX IN A SINGLE OSCILLATING 5-FUNCTION 

IMPURITY 

The scatterin g pr oblem of a single (5-function impurity with sinusoidal time dependence has been investigated by 
several authors [jjj, |43, |4Jj. We would like to summarize how to construct its Floquet scattering matrix in this 
Appendix. The system is described by the Hamiltonian 



where \i is the mass of the incident particle, while V s and Vd represent the strength of the static and the oscillating 
potential respectively. Using the Floquet formalism the solution of this Hamiltonian can be expressed as 
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H{x,t) 



^--r^ + [V s + V d cos(ojt + 4>)]S(x), 



(Al) 



OO 




= e 




(A2) 



n— — oo 



where Epi is the Floquet energy which take continuous values in the interval < Epi < huj. 
Since the potential is zero everywhere except at x = 0, ip n (x) is given by the following form 




x < 
x > 0, 



(A3) 
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where k„ 



VME- 



Fl 



nhuj)/h. The wave function ^ E Fi {x,t) is continuous at x = 0, 
A n + B n = C n + D n , 

and the derivative jumps by 



dx 



x=0+ 



dx 



x=0' 



2m 

= -j^[V s + V d coa(tJt + (f>)}y Em (0,t). 



(A4) 



(A5) 



Using Eq. (|A2|I this leads to the condition 



ik n (C n — D n A n + B n ) 
= ls (A n + B n ) + 7d (e-^A n+ i + e**An-i + e 
= ls {C n + D n ) + 7d (e-^C n+1 + e**C n -! + e~^D n+1 + e^D„-i), 



(A6) 



where 7 S = 2fiV s /h 2 and 7^ = fiVd/h 2 . After some algebra we have the following equation from Eqs. I|A4(1 and (|A6|) 

(3W-(i + r)-ir (/ + d-i \(a\ 

\c J \ (i + r)- 1 -(j + D^rJ [d J ' 1 j 

where 





yae^/ik-i js/ik-i w^/ik-! 
jae^/iko j s /ik W^/iko 



V 







(A8) 



/ 



and I is an infinite-dimensional square identity matrix. Eq. (| AT|) can also be expressed in the form |out) = M |in), 
where M connects the input coefficients to the output coefficients including the associated evanescent Floquet side- 
bands. In order to construct the scattering matrix we multiply an identity to both sides, K~ X K |out) = MK^ 1 K |in), 
where K nm = \fk~^8 nm . Then we have J ut_— MJ{ n , where J represents the amplitude of probability flux and 
M = KMK~ X . It should be mentioned that M is not unitary due to the evanescent modes included. If we keep only 
the propagating modes, we obtain the unitary scattering matrix S. 
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